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ABSTRACT 


A  new  type  of  acoustic  resonator  utilizes  alterations  of  the 
nonuniformity  to  achieve  different  resonance  frequencies.  Each  resonator 
is  designed  to  yield  frequencies  that  correspond  to  musical  notes.  The 
apparatus  are  remarkably  simple,  employing  piecewise  uniform  cross 
sectional  areas  that  can  easily  and  quickly  be  changed.  The  resonators 
are  thus  useful  as;  educational  demonstrations.  The  phenomenon  can  be 
understood  physically  and  quantified  perturbatively  by  the  energy  method 
of  Rayleigh,  or  equivalently  by  the  method  of  adiabatic  invariance. 
Improved  accuracy  is  obtained  by  matching  one-dimensional  solutions  of  the 
wave  equation  in  the  uniform  regions.  The  agreement  between  this  theory 
and  the  experimental  data  is  less  than  one  percent. 
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I.  INTRODUCTION 


In  this  chapter,  we  begin  with  the  background  of  and  motivation  for 
the  type  of  resonators  investigated  in  this  thesis.  Our  point  of  view  is 
that  the  resonators  are  primarily  useful  as  educational  demonstrations. 
We  discuss  these  demonstrations  in  detail. 

A.  BACKGROUND  AND  MOTIVATION 

The  modal  frequencies  of  any  resonator  depend  upon  the  geometry  of  the 
resonator  and  the  wave  speed  of  the  medium.  (If  the  medium  is  dispersive, 
then  the  wave  speed  depends  not  only  upon  the  medium  but  also  the 
frequency. )  Methods  of  changing  the  resonance  frequency  of  a  resonator 
fall  into  three  categories:  excitation  of  a  different  mode,  alteration  of 
the  wave  speed,  or  alteration  of  the  geometry.  A  dramatic  acoustics 
demonstration  of  the  first  is  the  “corrugahorn”  (Crawford,  1974),  and  of 
the  second  is  the  raising  or  lowering  of  the  characteristic  voice 
frequency  after  helium  or  sulfur  hexafluoride  is  inhaled.  Common 
demonstrations  of  the  third  involve  a  change  in  the  volume  (or  length)  of 
a  resonator.  It  is  not  well  known  that  another  means  of  changing  the 
geometry,  and  thus  the  frequency  of  a  mode,  is  to  introduce  or  alter  the 
nonun if ormitv  while  the  volume  is  held  constant.  In  the  case  of  the  one¬ 
dimensional  resonator,  the  nonuniformity  is  a  variation  in  the  cross 
sectional  area. 

We  have  constructed  simple  acoustic  resonators  whose  nonuniformity  can 
quickly  and  easily  be  altered  (Fig. I. A. 1).  The  sound  is  generated  by 
slapping  one  end  of  a  single-piece  resonator  by  hand  (Szeszol,  1987)  or 
with  a  rubber  stopper,  or  by  clapping  together  the  two  halves  of  a  double¬ 
piece  resonator  resting  on  a  smooth  surface.  The  fundamental  mode  of 
these  closed-open  resonators  is  excited.  By  slapping  the  other  end  of  a 
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single-piece  resonator,  or  changing  the  configuration  of  a  double-piece 
resonator,  we  effectively  alter  the  nonuniformity  and  thus  the  frequency. 
The  apparatus  are  useful  as  an  educational  demonstration  of  the 
effects  of  a  nonuniformity  on  standing  waves.  Each  resonator  is  designed 
such  that  the  differences  in  frecpiency  correspond  to  musical  intervals. 
This  aides  in  the  listeners'  recognition  of  the  change  in  frequency,  makes 
the  effect  more  dramatic,  and  indicates  the  utility  of  theory.  In  the 
next  section  we  describe  the  demonstrations  in  detail.  Although  the 
phenomenon  can  be  appreciated  on  an  elementary  level  as  a  result  of 
symmetry-breaking,  the  explanation  of  the  effect  is  not  elementary.  In 
fact,  our  experience  is  that  no'-.^cousticians  and  even  many  acousticians 
are  unable  to  explain  the  effect.  In  Chapter  II,  we  give  three  different 
theoretical  treatments  of  the  problem:  matching  solutions  of  the  wave 
equation,  applying  the  energy  method  of  Rayleigh,  and  applying  the  theorem 
of  adiabatic  invariance.  The  latter  two,  which  are  perturbative,  offer 
physical  explanations  for  the  effect.  In  Chapter  III,  we  discuss  the 
design  of  the  resonators,  and  in  Chapter  IV,  quantitatively  compare  the 
measured  frequencies  to  the  predicted  values. 
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Pig. I. A. 1  Variably  nonuniform  acoustic  resonators:  (a)  single-piece 
"slapper",  and  (b)  double-piece  "clapper". 


B.  DEMONSTRATIONS 


The  following  is  an  effective  way  to  conduct  the  demonstrations. 
First,  a  single-piece  resonator  with  uniform  inner  diameter  is  slapped  by 
hand  on  one  end.  The  sound  can  be  easily  heard  in  a  large  lecture  hall. 
The  other  end  is  then  slapped  to  show  that  the  same  frequency  is  produced, 
as  required  by  symmetry.  The  symmetry  is  broken  in  the  next  resonator, 
which  has  a  discontinuity  in  dicimeter  (Fig.  I. A.  la) .  The  two  frequencies 
are  now  different.  Three  such  resonators  are  demonstrated,  each  with  a 
successively  greater  relative  difference  in  diameter.  The  first  resonator 
is  designed  to  have  a  frequency  difference  corresponding  to  a  musical 
interval  of  a  major  third,  the  second  resonator  a  perfect  fifth,  and  the 
third  resonator  an  octave. 

How  can  an  additional  frequency  be  obtained  from  a  resonator?  If  a 
two-pitch  resonator  is  imagined  to  be  cut  in  half  along  the  axis 
(Fig. I .A. lb) ,  there  are  now  three  distinct  closed-open  configurations,  and 
thus  three  different  frequencies  will  occur.  Three  such  resonators,  which 
correspond  to  halving  the  above  three,  are  then  demonstrated.  Rather  than 
slapping  the  ends,  it  is  now  convenient  to  clap  together  the  two  halves 
while  one  set  of  ends  rests  on  a  smooth  surface.  The  first  resonator 
yields  the  notes  1,2,3  of  a  major  scale,  the  second  resonator  1,3,5,  and 
the  third  resonator  1,5,8  (in  which  the  final  frequency  is  an  octave  above 
the  first ) . 

Of  the  four  possible  closed-open  configurations  above,  two  are 
degenerate.  To  obtain  a  four-pitch  resonator,  we  construct  the  two  halves 
to  be  dissimilar.  The  symmetry  is  then  completely  broken.  The  final 
demonstration  involves  two  four-pitch  resonators  that  are  designed  to 
yield  the  eight  notes  of  a  musical  scale  (in  which  the  first  and  final 
notes  differ  by  an  octave) . 
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In  all  of  the  above  demonstrations,  each  pitch  is  produced  three  to 
five  times  in  rapid  succession,  so  that  the  frequency  can  be  implanted  in 
the  listeners'  minds.  This  facilitates  the  comparison  with  successive 
frequencies . 
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II.  THEORY 


In  this  chapter,  we  develop  three  different  approaches  to  determine 
how  the  frequencies  of  a  resonator  depend  upon  nonuniformity.  Several 
special  cases  are  investigated  in  detail. 

A.  STANDARD  APPROACH 

We  consider  the  acoustic  resonator  shown  in  Fig. II. A. 1.  The  problem 
is  to  determine  the  frequencies  of  the  "quasilongitudinal"  modes.  The 
term  refers  to  the  fact  “-hat  the  flow  for  these  modes  is  not  purely 
I'-ngitudinal ,  due  to  the  discontinuity  in  diameter.  In  this  chapter,  we 
assume  that  the  length  of  the  resonator  is  much  greater  than  the 
characteristic  diameter,  so  that  the  end  correction  can  be  neglected.  The 
solution  of  the  problem  when  Aj  =  A,,  (or,  equivalently,  a  =  0  or  a  =  1) 
is  elementary.  The  angular  frequency  of  the  n''’  longitudinal  mode  is 

'  (II. A. 1) 

where  c  is  the  speed  of  sound,  and  where  the  wavenumber  is 

Jc.  -  -i  (2„-l)  ,  „-l,2,3 .  (II.ft.2) 


To  solve  the  problem  when  A2  A,,  we  will  assume  that  the  flow  is 
purely  longitudinal  except  in  a  narrow  region  about  the  step.  The  problem 
then  reduces  to  boundary  condition  matching  of  the  svjlutions  of  the  one¬ 
dimensional  wave  equation,  as  we  now  show. 

We  begin  with  the  linearized  one-dimensional  ideal  fluid  equations. 
The  linearized  Euler's  equation,  which  is  a  statement  of  Newton's  Second 
Law,  is 
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(II. A. 3) 


P 


8v 

‘>at 


where  is  the  equilibrium  density,  v  is  the  velocity  and  5p  is  the 
deviation  of  the  pressure  from  the  equilibrium.  The  linearized  equation 
of  continuity,  which  is  a  statement  of  mass  conservation,  is 

(II.A.4, 


where  6p  is  the  deviation  of  the  density  from  equilibrium.  Ecjuations 
(II. A. 3)  and  (II. A. 4)  are  valid  if  v,  6p,  and  6p  are  small.  Under  these 
conditions  the  pressure  and  density  are  related  thermodynamically  by 


5p  -  c^6p  . 


(II. A. 5) 


By  eliminating  5p  and  v  in  (II. A. 3),  (II. A. 4),  and  (II. A. 5),  one  can 
readily  derive  the  wave  equation  for  5p: 


-  c^  — «p  -  0  . 

dt^  dx^ 


(II. A. 6) 


The  velocity  v  satisfies  the  same  equation.  The  wave  equation  (II. A. 6)  is 
valid  in  region  1  (to  the  left  of  the  step  in  Fig.  II. A. 1)  and  in  region 
2  ( to  the  right ) . 

To  consider  standing  waves,  we  set  5p  =  p' (x)exp(iu)t) ,  where  w  is  as 
yet  unknown.  The  wave  equation  (II. A. 6)  then  becomes 
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The  boundary  conditions  at  x  =  oL  are  that  the  pressure  and  volume 
velocity  Av'  are  both  continuous,  where  A  is  the  cross-sectional  area.  By 
(II.  A.  8),  the  latter  is  equivalent  to  the  continuity  of  hdp' /dx.  Note 
that  dp' /dx  is  thus  discontinuous  for  our  system.  By  imposing  the 
boundary  conditions  on  (II. A. 9)  and  (II. A. 10),  and  dividing  the  resultant 
equations  by  each  other  in  order  to  eliminate  the  amplitudes  a  and  b,  we 
find  the  wavenumber  condition 
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(II. A, 11) 


tan  (kLa)  tan  [AIi(l  -  a)  ] 


A, 


The  problem  is  to  determine  the  values  of  k  that  satisfy  this  equation. 
The  frequency  w  is  then  given  simply  by  the  dispersion  law  w  =  ck. 

For  the  special  case  Aj  =  A,  (no  step),  (II. A. 11)  becomes 

tan  (kLa)  tan  (kL(l  -  a)  ]  -  1  .  (II. A. 12) 

This  should  imply  the  elementary  solution  (II. A. 2).  It  is  not  even 
mathematically  obvious  that  the  k  solution  of  (II. A. 12)  is  independent  of 
a,  although  we  know  physically  that  it  must  be.  To  investigate  this 
special  case,  it  is  convenient  to  redefine  a  by  letting  a  -♦  a  +  1/2. 
Equation  (I I. A. 12)  then  becomes 

tan  (kL(l/2  +  a)  ]  tan  (kI,(l/2-a))-l  .  (II. A.  13) 


Using  the  identity 


tan(a  0) 


tan  a  +  tan  0 
1  -  tan  a  tan  0  ' 


(II. A. 14) 


and  simplifying,  we  find  tan^(kl,/2)  =  1.  This  implies  k  =  (7r/2L)  (2n-l) , 
n  =  1,2,3,...,  which  indeed  is  (II. A. 2). 

Equation  (II. A. 11)  is  invertible  if  a  =  1/2,  i.e.,  when  the  step  is 
midway  between  the  ends  of  the  resonator.  The  fractional  frequency 
of  the  fundamental  relative  to  its  unperturbed  (Aj  =  A,)  value  is  then 
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(II. A. 15) 


In  Fig. II. A, 2,  the  fractional  frequency  shift  {w-Wo)/Wo  is  plotted  as  a 
function  of  the  relative  area  difference  (Aj-A, )  /  (Aj+A, )  . 
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Pig. 11. A. 1  Geometry  of  the  acoustic  resonator.  The  cross  sectional  area 
equals  A,  for  OsxSaL,  and  for  aL<x<L,  where  0<a<l. 
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B.  SPECIAL  CASES  OF  THE  STANDARD  APPROACH 


In  this  section,  we  consider  three  special  cases  of  the  condition 
(II. A. 11)  obtained  by  the  standard  approach.  The  first  case  is  a  small 
step,  the  second  is  a  large  narrowing  step  (Aj  «  A,  in  Fig.II.A.l),  and 
the  third  is  a  large  widening  step  (A^  »  A,).  These  cases  are  interesting 
for  a  variety  of  reasons,  the  moat  important  being  that  they  can  be 
understood  qualitatively.  The  cases  will  be  examined  independently  and 
thoroughly  in  the  following  paragraphs. 

We  now  consider  a  small  step,  or  the  case  where  A2  is  approximately 
equal  to  A,.  More  precisely,  this  is 


1  <  <  1 


(II.B.l) 


We  refer  to  this  as  Ravleioh^a  case,  because  he  derived  a  formula  for 
frequency  shifts  due  to  slight  nonuniformities  (Section  II. C). 

The  k  solution  of  the  condition  (II. A. 11)  will  be  approximately  equal 
to  the  elementary  (Aj  =  A,)  case  (II.A.2).  Hence,  we  set 

k  -  k„(l  +  c)  ,  (II. B. 2) 


where  |el  «  1.  Note  that 


(II. B. 3) 


where  6k  is  defined  as  k  -  k„.  We  wish  to  determine  e.  Once  this  is 
found,  k  is  determined  from  (II.B.3),  and  w  is  determined  by 
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(II. B. 4) 


6<o  _  bk 

”  K 


where  Su  =  u  -  a)„.  The  relationship  (II.B.4)  follows  exactly  from  the 
dispersion  law  u  =  ck.  The  expression  (II. B. 2)  is  substituted  into  the 
condition  (II. A. 11),  which  is  then  expanded  to  the  first  order  in  e.  The 
expansion  of  the  tangents  is  given  by 

tan  [(0(l  +  e))]  -  tan0  +  e0  sec^0  ,  (II. B. 5) 

which  is  valid  for  jej  <<  1.  After  employing  some  elementary 

trigonometric  identities,  and  simplifying,  we  find 

.  sin[ax(2n-l)l  (II. 8.6) 

w„  jt(2n-l)  Aj 


where  we  have  used  (II. B. 3)  and  (II.B.4).  For  the  fundamental  (n  =  1), 
this  is 


3<ii  _  sin(ait)  (4^-1) 

0)i  Ji  A, 


(II. B. 7) 


The  frequency  shift  fiw  vanishes  when  a  =  0  or  a  =  1,  and  is  an  extremum 
when  a  =  1/2.  We  will  derive  (II. B. 5)  physically  in  later  sections  (II. C 
and  II .D) . 

Next,  we  specialize  the  condition  (II. A. 11)  to  the  case  of  a  large 
narrowing  step,  Aj  «  A,.  The  system  then  becomes  a  Helmholtz  resonator. 
For  the  fundamental,  we  can  use  the  approximation  tan  $=  6  in  the 
condition  (II. A. 11),  resulting  in 

k^L^a{l-a)  -  ^  .  (II. B. 8) 
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which  is  valid  for  Aj  «  .  By  using  the  dispersion  law  o)  =  ck,  and 

solving  for  the  frequency,  we  find 


u 


2 


1 

L^a{l-a) 


(II. B. 9) 


We  will  show  in  Section  II. E  that  this  is  equivalent  to  the  well-known 
result  obtained  by  directly  deriving  the  frequency  of  a  Helmholtz 
resonator. 

Finally,  we  consider  a  large  widening  step,  Aj  »  A,.  In  this  case, 
the  condition  (II. A. 11)  becomes 

[tan  kLa]  [tankL(l-a)]  -  «  .  (II.B.IO) 


This  implies 


k 


n 

2  La 


(2n-l) 


n-  1,2,3,  ..  . 


(II.B.ll) 


or 


(II.B.12) 


If  a  =  1/2,  then  (II.B.ll)  and  (II.B.12)  become  identical.  The  density  of 
modes  is  thus  half  as  much  as  when  a  1/2.  We  will  show  in  Section  II. F 
that  this  degeneracy  is  lifted  when  A^/A,  ^ 
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C.  RAYLEIGH'S  APPROACH 

In  this  section  we  determine  the  frequency  shift  of  a  mode  in  a 
closed-open  resonator  with  a  small  arbitrary  nonuniformity  (Fig. II.C. 1) . 
The  calculation  is  similar  to  that  of  Rayleigh  (1896).  The  results  can  be 
easily  transformed  to  the  cases  of  closed-closed  and  open-open  resonators. 

The  instantaneous  potential  and  kinetic  energy  densities  of  the 
unperturbed  mode  will  be  denoted  by  pe(x,t)  and  ke(x,t),  respectively. 
The  nonuniformity  6A(x)  causes  the  total  instantaneous  potential  energy  PE 
to  change  by  the  amount 

6PE  -  ipeSA  *  A5pe)dx  ,  (II.C.l) 


where  6pe  is  to  be  determined.  The  change  in  the  total  instantaneous 
kinetic  energy  KE  is  similar,  except  that  there  is  an  additional  term  due 
to  a  possible  frequency  shift: 

5KE  -  f‘  (ke5A  +  ASke)dx  +  2—KE  .  (II.C. 2) 

Jo  oi 

The  factor  of  two  arises  because  KE  is  proportional  to  the  square  of  w. 
The  potential  energy  density  is  proportional  to  the  square  of  the  density 
amplitude  p’ .  Furthermore,  the  fractional  change  in  density  amplitude 
equals  the  negative  of  the  fractional  change  in  volume,  by  continuity  of 
mass.  Hence, 

-  iIEI  -  -2—  .  (II.C. 3) 

pe  p'  A 


The  fractional  change  in  kinetic  energy  density  is  similarly 
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5ke  _  2  _  _2 

ke  IT  ~A 


(II. C. 4) 


where  v  is  the  velocity  amplitude.  Substituting  (II. C. 3)  and  (II. C. 4) 
into  (II.C.l)  and  (II. C. 2),  respectively,  gives 

SPE  -  - pe  SA  dx  ,  (II.C.5) 


6KE  -  ke  6A  dx  *  2—  KE  .  (II. C. 6) 

Jo  w 

For  a  mode  executing  simple  harmonic  motion  in  which  all  mass  elements  are 
in  phase,  the  maximum  total  potential  energy  must  equal  the  maximum  total 
kinetic  energy.  We  thus  set  =  6KE^^,  and  solve  for  the  frequency 

shift.  The  result  is 

i  r  '  (II. c. 7) 

where  E  is  the  total  energy  of  the  unperturbed  mode.  Finally,  we 
substitute  the  explicit  forms  of  the  maximum  energy  densities: 

pe„„  -  ^  cos^kx  .  (II. C. 8) 


where  the  wavenumber  of  the  mode  is  k  =  (2n-l)7r/2L,  n  =  1,2,3...,  for  a 
closed-open  resonator.  The  final  result  is 
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6a) 

"aT 


(II.C.IO) 


^  coa(2/cx)  6A  dx  , 

AL  Jo 

which  is  equivalent  to  the  result  of  Rayleigh  (Rayleigh,  1945). 

The  effect  can  be  understood  qualitatively.  Suppose  that  a  slight 
constriction  of  small  length  exists  at  some  location  along  the  resonator. 
If  the  location  is  not  a  pressure  node  of  the  mode  of  interest,  the 
potential  energy  of  the  mode  will  be  altered.  That  the  potential  energy 
density  exists  over  a  smaller  volume  causes  the  potential  energy  to 
decrease.  However,  the  density  amplitude  of  the  mode  increases  as  a 
result  of  mass  conservation,  causing  an  increase  in  the  potential  energy. 
The  fractional  decrease  in  volume  is  identical  to  the  fractional  increase 
in  density.  Because  the  potential  energy  density  is  proportional  to  the 
scmare  of  the  density  amplitude,  the  magnitude  of  the  increase  in  energy 
is  twice  the  magnitude  of  the  decrease.  The  net  effect  is  thus  an 
increase  in  potential  energy.  A  constriction  at  a  location  that  is  not  a 
pressure  node  therefore  causes  the  stiffness  of  the  mode  to  increase, 
which  will  raise  the  frequency.  The  opposite  occurs  if  there  is  a  bulge 
rather  than  a  constriction.  There  is  no  effect  on  the  stiffness  if  an 
area  change  is  located  at  a  pressure  node.  We  can  check  the  result  by 
considering  a  constriction  at  a  closed  end  of  a  resonator.  The  frequency 
indeed  increases  in  this  case,  because  the  constriction  is  equivalent  to 
a  shortening  of  the  length  of  the  resonator. 

To  consider  the  effect  upon  the  kinetic  energy,  we  first  note  that  if 
the  location  of  an  area  change  is  at  a  velocity  node,  then  the  kinetic 
energy  is  unaltered.  Otherwise,  there  is  a  competition  between  volume  and 
amplitude  changes  analogous  to  the  potential  energy  case.  For  a 
constriction,  the  velocity  increases  due  to  mass  conservation,  and  the 
magnitude  of  the  increase  in  kinetic  energy  is  twice  the  magnitude  of  the 
decrease  due  to  the  volume  decrease.  The  net  effect  is  that  a 
constriction  at  a  location  that  is  not  a  velocity  node  causes  the  kinetic 
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energy  to  increase.  There  is  thus  an  increase  in  the  inertia  of  the  mode, 
and  the  frequency  decreases. 

In  conclusion,  a  small  constriction  at  a  velocity  node  increases  the 
stiffness  and  has  no  effect  upon  the  inertia,  so  the  frequency  increases. 
At  a  pressure  node,  a  small  constriction  increases  the  inertia  and  has  no 
effect  upon  the  stiffness,  so  the  frequency  decreases.  The  effects 
compete  if  the  area  alteration  is  at  neither  a  velocity  nor  pressure  node. 
In  the  case  where  the  cross  sectional  area  continuously  changes  along  the 
length  of  the  resonator  such  that  the  deviation  from  the  mean  value  is 
always  small,  the  frequency  shift  is  the  integral  over  differentiax  shifts 
due  to  area  changes  of  differential  length. 

In  addition  to  being  physical,  Rayleigh's  energy  method  is  in  general 
very  useful  as  a  means  of  determining  the  frequency  shifts  due  to  small 
nonuniformities.  For  example,  it  has  been  applied  to  the  modes  of  a  bar 
whose  ends  have  additional  mass  due  to  transduction  coils  (Garrett,  1990). 

The  Rayleigh  formula  (II.C.IO)  exhibits  a  symmetry  that  is  important 
in  the  design  of  our  resonators.  If  the  closed  and  open  ends  are 
switched,  note  that  6a)/oj  simply  changes  sign.  There  thus  exists  a 
reflection  antisymmetry.  If  a  resonator  is  composed  of  two  halves,  there 
are  then  four  possible  closed-open  orientations  and  thus  four  frequencies 
(which  may  not  be  distinct).  As  a  result  of  the  reflection  antisymmetry, 
however,  the  relationships  of  the  frequencies  is  not  arbitrary. 

To  quantify  this,  we  define  the  functional 

7,  -  cos(2/tx)  5A,  dx  ,  (II.C.ll) 

and  similarly  72,  where  5A,(x)  and  6A2(x)  are  the  variations  in  cross 
sectional  area  of  the  halves  that  compose  the  resonator.  The  total 
variation  is  5A  =  5A,  +  6A2.  The  fractional  frequency  shifts  Sui/u 
corresponding  to  the  four  orientations  are  then  *(71+72)  and  ±(7,-72). 
Without  loss  of  generality,  we  assume  that  7,, 72  2  0  and  72  2  7,.  For  a 
fixed  function  5A2(x),  and  thus  a  fixed  value  of  y^,  we  consider  different 
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functions  6A,  (x)  such  that  0  <  7,  s  72.  This  explores  all  possible 
relationships  cimong  the  four  f recjuencies.  The  results  are  shown  in 
Fig. II. C. 2.  The  orientations  of  the  two  halves  are  labeled  by  arrows.  We 
have  arbitrarily  labeled  the  lowest  frequency  orientation  by  two  down 
arrows.  The  degenerate  cases  in  Fig. II. C. 2  can  occur  as  follows.  If  one 
of  the  halves  is  uniform,  then  5A,  =  0  and  thus  7,  =  0.  There  are  then 
only  two  frequencies.  If  both  halves  are  identical,  then  5A,  =  6A2  and 
thus  7,  =72-  There  are  then  only  three  frequencies.  These  special  cases 
assume  that  72  ^  0. 

Finally,  we  specialize  the  Rayleigh  formula  (II.C.IO)  to  the  case  of 
the  fundamental  mode  of  our  piecewise  uniform  resonator  (Fig. II. A. 1) . 
After  integrating  and  simplifying,  we  find 

_6w  _  sin  TO  _  (II.C.12) 

(ti  A  n  ' 

This  is  valid  if  the  relative  difference  between  A,  and  Aj  is  small.  The 
reference  area  A  can  be  chosen  as  any  value  that  roughly  equals  A,  (or  Aj). 
(Deviations  in  (II.C.IO)  or  (II.C.12)  due  to  the  choice  of  A  are  beyond 
the  accuracy  of  the  result.)  Comparison  of  (II.C.12)  witn  the  special 
case  (II.B.7)  of  the  standard  result  shows  that  the  two  are  equivalent  to 
the  order  of  the  approximation.  For  a  =  1/2  and  A  =  (Ai+A2)/2,  a  plot  the 
Rayleigh  result  is  shown  in  Fig. II. A. 2.  Note  that  it  closely  approximates 
the  standard  result  over  a  broad  range  of  differences  in  the  areas. 
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Fig.II.C.l  Closed-open  resonator  with  a  small  arbitrary 
nonuniformity. 


Fig. II. C. 2.  Fractional  frequency  splittings  of  an  open-closed 
resonator  composed  of  two  halves. 


D.  ADIABATIC  INVARIANCE  APPROACH 

In  the  expression  (II. C. 7),  we  observe  that  the  frequency  shift  due 
to  a  perturbation  in  cross  sectional  area  at  a  location  depends  only  upon 
the  relative  difference  of  the  potential  and  kinetic  energy  densities  at 
that  location.  One  might  wonder  if  this  is  a  general  result  which  applies 
in  any  geometry  (including  higher  dimensions).  By  employing  adiabatic 
invariance,  we  now  show  that  this  is  indeed  true. 

We  consider  any  mode  of  an  arbitrary  cavity.  To  neglect  end  effects 
due  to  possible  openings,  we  assume  that  either  the  cavity  is  completely 
closed  or  that  periodic  boundary  conditions  occur.  (The  latter  is 
appropriate  for  the  case  of  our  one-dimensional  closed-open  resonator.) 
Suppose  that  the  mode  has  energy  E  and  frequency  w,  and  that  dissipation 
is  neglected.  We  now  imagine  that  an  external  agent  slowly  deforms  the 
surface  of  the  cavity  by  a  small  amount,  such  that  the  displacement  at 
every  point  is  perpendicular  to  the  surface.  By  the  theorem  of  adiabatic 
invariance,  energy  remains  in  a  single  mode.  Because  the  motion  is  simple 
harmonic,  the  adiabatic  invariant  is  the  quotient  of  the  energy  and  the 
frecjuency  (Crawford,  1990).  Hence,  in  the  perturbative  case. 


_  6E 

~zir  ~E' 


(II.D.l) 


where  w  +  6w  is  the  new  frequency  and  E  +  6E  is  the  new  energy. 

The  change  in  energy  5E  equals  the  work  6W  done  by  the  external  agent. 
If  p  is  the  time-averaged  pressure,  dS  is  an  element  of  surface  area,  and 
is  the  displacement,  the  work  is 
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5W  -  ^pdSSi 


(II. D. 2) 


where  the  integration  ia  over  the  surface  of  the  cavity.  To  second  order 
in  the  amplitude  of  the  sound  field,  the  time-averaged  pressure  is  (Landau 
and  Lifshitz,  1959) 

P  ■  Po  ^  +  const  ,  (II. D. 3) 

where  the  radiation  pressure  on  the  surface  of  the  cavity  is 

P«rf  -  (PCm.,  -  )  •  (II. D. 4) 

In  these  expressions,  p,  is  the  equilibrium  pressure,  and  pe„„  and  ke„„  are 
the  maximum  potential  and  kinetic  energy  densities  at  a  point.  The  second 
term  in  the  radiation  pressure  (II.D.4)  is  the  familiar  Bernoulli  pressure 
that  also  arises  in  incompressible  flow.  The  first  term  is  an  additional 
pressure  that  arises  due  to  the  compressibility.  In  the  time-averaged 
pressure  (II.D.3),  the  constant  is  necessary  in  order  to  keep  the  total 
mass  fixed,  because  there  is  also  a  second  order  time-averaged  density 
which  is  in  general  nonzero. 

The  part  of  the  work  (II.D.2)  due  to  the  pressure  p„  +  const 
corresponds  to  a  change  in  the  overall  density,  and  therefore  to  a  change 
in  temperature.  We  ignore  this  part  of  the  work  because  we  desire  the 
frequency  shift  due  to  geometrical  and  not  temperature  alteration.  By 
combining  (II.D.l),  (II.D.2),  (II.D.3),  and  (II.D.4),  we  find  that  this 
frequency  shift  is 


24 


6a; 

"w* 


f  (pe 

m«x  -  )  ds  . 


(II. D. 5) 


To  apply  this  to  the  case  of  a  one-dimensional  resonator  that  is  initially 
uniform,  we  set  dS5$  =  dx£A,  where  dx  is  an  element  of  length  of  the 
resonator  and  6A(x)  is  the  variation  in  the  cross  sectional  area.  The 
result  (II. D. 5)  is  then  identical  to  the  intermediate  result  (II. C. 7)  in 
Rayleigh's  case. 

Adiabatic  invariance  offers  an  elegant  explanation  of  the  frequency 
shift  due  to  nonuniformity  of  a  resonator.  In  fact,  we  can  view  the 
frequency  shift  as  an  indirect  demonstration  of  adiabatic  invariance.  On 
a  practical  level,  the  method  predicts  the  frequency  shift  when  an 
arbitrarily  complicated  geometry  is  altered  a  small  amount  at  a  point. 
All  that  is  required  is  the  relative  difference  of  the  potential  and 
kinetic  energy  densities  at  that  point. 
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B .  HELMHOLTZ  RESONATOR 


In  the  limit  Aj  «  A,,  our  piecewise  uniform  resonator  (Fig. II. A. 1) 
becomes  a  Helmholtz  resonator.  In  this  section,  we  directly  derive  the 
frequency  of  this  resonator,  and  show  that  the  result  is  equivalent  to 
the  limiting  case  of  our  general  result. 

The  geometry  of  a  Helmholtz  resonator  is  shown  in  Fig.II.E.l.  The 
frequency  can  be  determined  by  applying  Newton's  Second  Law  to  the  mass  of 
the  fluid  in  the  neck.  The  amount  of  this  mass  is 

m  -  p  h  A  ,  (lI.E.l) 


where  p  is  the  density.  If  the  volume  hA  of  the  neck  is  much  less  than 
the  volume  V  of  the  cavity,  the  fluid  in  the  neck  then  oscillates  rigidly 
and  the  fluid  in  the  cavity  contracts  and  expands  uniformly.  The  excess 
force  exerted  by  the  fluid  in  the  cavity  on  the  fluid  in  the  tube  is 

F  -  A  p'  ,  (II. E. 2) 


where  p'  is  the  excess  pressure  in  the  cavity.  Using  p'  =  c^p',  where  p' 
is  the  excess  density  in  the  cavity,  and 


p'  _  _  Ax 

~  ~  V  ~  V 


(II. E. 3) 


we  can  express  (II.E.2)  as 


F  - 


A^x 

V 


(II. E. 4) 


Substituting  this  and  (lI.E.l)  into  F  =  md^x/d^t  gives 
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(II. E. 5) 


d*x  c^A 

dP 


0  , 


from  which  we  find  the  frec[uency; 


u 


2 


c^A 


(II. E. 6) 


This  is  the  standard  expression  for  the  frequency  of  a  Helmholtz 
resonator.  To  apply  (II.E.6)  to  the  geometry  of  our  piecewise  uniform 
resonator  (Fig. II . A. 1 ) ,  we  set  A  =  Aj,  h  =  (l-a)L,  and  V  =  A,aL.  These 
yield 


O) 


2 


1 

L^a  ( 1  -  a ) 


(II. E. 7) 


This  is  identical  to  (II. B. 9),  which  is  our  general  result  evaluated  in 
the  limit  Aj  «  A,. 

We  can  express  (II. E. 7)  in  dimensionless  variables  by  employing  the 
unperturbed  frequency,  which  is  w,  =  »rc/2L.  Then, 


<i)  _  2  1 

“o  ( 1  -  a ) 


(II. E. 8) 


This  is  plotted  in  Fig. II. A. 2  together  with  the  standard  result  and  the 
Rayleigh  approximation.  The  standard  result  connects  the  Helmholtz  and 
Rayleigh  cases. 
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Pig.II.E.l  Geometry  of  a  Helmholtz  resonator. 
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P.  WEAKLY  COUPLED  imiFORM  RESONATORS 

In  this  section  we  consider  the  case  of  a  large  widening  step, 
kj  »  A,  in  Fig. II. A. 1.  We  found  in  (II.B.IO)  and  (II.B.ll)  that  when 
Aj/A,  =  00  in  the  condition  (II. A. 11),  the  frequency  spectrum  is  the  union 
of  the  two  simple  spectra  of  uniform  resonators  of  length  aL  and  (l-a)L. 
In  this  limiting  situation,  therefore,  we  effectively  have  two  uncoupled 
uniform  resonators.  This  is  reasonable  physically;  the  resonator  with  the 
small  cross  sectional  area  A,,  which  has  one  closed  end,  has  its  other  end 
effectively  completely  open.  Similarly,  the  resonator  with  the  large 
cross  sectional  area  A^,  which  has  one  open  end,  has  its  other  end 
effectively  completely  closed.  Fig.II.F.l  shows  the  frequency  spectrum  as 
a  function  of  the  step  location  a,  where  0  <  a  <  1.  The  dashed  curves 
correspond  to  the  limiting  case  Aj/A,  =  in  which  there  are  two  uncoupled 
uniform  resonators. 

When  1  «  Aj/A,  <  <»,  these  uniform  resonators  become  weakly  coupled. 
In  this  case,  we  determine  the  spectrum  by  numerically  solving  (II. A. 11). 
The  solid  curves  in  Fig.II.F.l  correspond  to  Aj/A,  =  20.  The  phenomenon 
of  mode  level  repulsion  occurs  in  the  regions  where  the  uncoupled  (dashed) 
curves  intersect:  the  weak  coupling  lifts  the  degeneracy  here.  The 
situation  is  essentially  the  same  as  in  the  elementary  case  of  two  weakly 
coupled  identical  oscillators.  The  piecewise  uniform  resonator  in  fact 
offers  one  of  the  simplest  examples  of  mode  level  repulsion  in  a 
continuum. 

Finally,  we  comment  on  the  nature  of  the  spectrum  in  Fig.II.F.l  when 
Aj  is  not  much  greater  than  A,.  When  Aj  =  A,,  the  solid  curves  become 
horizontal  lines  at  the  values  w/w,  =1,3,5,...  .  (This  is  the  elementary 
case  of  a  uniform  resonator.)  When  A^/A,  «  1,  we  again  have  two  weakly 
coupled  resonators,  but  now  the  boundary  conditions  for  each  resonator  are 
not  mixed.  (One  resonator  is  approximately  open-open;  the  other  closed- 
closed.)  The  uncoupled  spectrum  curves  are  thus  hyperbolas  as  before,  but 
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now  they  have  the  values  w/w,  =  2,4,6,...  when  a  =  0  or  a  =  1.  The 
fundamental  is  the  Helmholtz  mode,  and  the  higher  modes  will  exhibit  the 
level  repulsion  effect  mentioned  above.  Note  that  the  curves  must  still 
have  the  values  u/w,  =  1,3,5,...  when  a  =  0  or  a  =  1,  however.  Near  these 
values  of  a  the  coupling  ceases  to  be  weak;  the  boundary  conditions  are  no 
longer  approximately  unmixed. 

The  numerical  solution  that  yields  the  solid  curves  in  Fig.II.F.l  is 
not  standard  and  should  thus  be  explained.  The  problem  can  be  posed  as 
the  determination  of  values  of  such  that  f  =  0,  where 


-tan{JLJLa)  tan[”J£-(l  -  a)] 
o).  2  w.  2 


A, 


There  are  an  infinite  number  of  such  zeros  (or  roots),  where  the  first 

lies  between  (jii/oi,  =  0  and  =  2,  the  second  between  2  and  4,  the  third 

between  4  and  6,  etc.  The  bisection  method  is  thus  preferable  to  Newton's 
method.  A  problem  with  the  bisection  method,  however,  is  that  f  has 
poles.  To  the  bisection  method,  these  appear  as  roots.  The  problem 
cannot  be  removed  by  redefining  f  such  that  the  tangents  appear  in  the 
denominator,  because  now  the  zeros  of  the  tangents  give  rise  to  poles. 
Our  trick  is  to  multiply  f  by  functions  that  have  simple  zeros  where  f  has 
its  (simple)  poles.  Specifically,  to  determine  the  five  curves  in 
Fig.II.F.l,  we  use  bisection  to  solve  PiPjf  =  0,  where  the  prefactors  are 

p,  -  (  JLa  -  1)  (  J£_a  -  3)...(_iLa  -  9)  , 

Pj  -  [  J£_(l  -  a)  -  1)  (  JfL(l  -  a)  -  3)..  .  [  J!L(1  -  a)  -  9)  . 

Wo  W„ 

Note  that  this  technique  does  not  alter  the  values  of  the  roots. 


30 


Frequency 


Fig.II.F.l  Frequency  spectrum  of  a  piecewise  uniform  resonator  when 
Aj/A,  =  20  (solid  curves),  and  A^  A,  =  “  (dashed  curves). 
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III.  DESIGN  AND  CONSTRUCTION 


In  this  chapter,  we  discuss  the  design  and  construction  of  the 
resonators  such  that  the  ratios  of  the  frequencies  will  approximately 
equal  musical  values.  Despite  the  simple  appearance  of  the  resonators,  we 
will  show  that  neither  the  design  nor  construction  are  easily 
accomplished. 

A.  DESIGN  OF  THE  2-PITCH  RESONATORS 

By  a  "two-pitch"  resonator  we  mean  a  single-piece  resonator  that  is 
slapped  by  hand  at  one  end  and  then  the  other.  Because  this  switches  the 
boundary  conditions  on  the  acoustic  field,  we  have  seen  that  a  nonuniform 
resonator  will  in  general  have  two  different  frequencies.  In  this  section 
we  discuss  the  design  of  these  resonators  (called  "slappers")  such  the 
ratio  of  the  frequencies  corresponds  to  a  musical  interval. 

An  octave  is  an  interval  from  one  frequency  to  twice  that  frequency. 
In  the  equal-tempered  musical  scale,  an  octave  is  split  into  twelve 
intervals  that  are  equally  spaced  with  respect  to  the  logarithm  of  the 
frequency.  This  is  necessary  if  a  sequence  of  frequencies  (or  notes)  is 
to  be  perceived  as  essentially  the  same  regardless  of  the  key  (or  first 
frequency  of  the  scale).  Let  f„,  n  =  0, 1,2, . . . , 12,  be  the  frequencies  of 
a  chromatic  scale.  Equal  temperment  implies 

log  f„  -  log  ^  (log  ^12  -  log  f„)  .  (III.A.l) 

Substituting  f,2  =  2f„,  simplifying,  and  exponentiating,  yields 

-  _  2^  ,  (III. A. 2) 
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The  relative  frequencies  fn/fo  tabulated  in  Table  III.A.l.  We  have 
labeled  the  frequencies  not  by  n,  but  by  the  notes  corresponding  to  the 
major  scale,  which  is  the  sequence  1, 2 , 3, 4, 5 , 6, 7 , 8.  We  will  design 
slappers  whose  frequency  ratios  are  3,  5,  and  8.  The  two  generated  notes 
can  thus  be  denoted  as  1-3,  1-5,  and  1-8,  which  is  how  we  will  label  the 
resonators . 

The  simplest  nonuniformity  to  impose  is  a  step  in  an  otherwise  uniform 
resonator  (Fig. II. A. 1)  We  investigated  this  case  in  Section  II. A.  If  the 
end  correction  is  neglected,  the  relationship  between  frequency  and  the 
geometry  is  given  by  Equation  (II. A. 11).  For  a  fixed  discontinuity  in  the 
cross  sectional  area,  the  greatest  frequency  shift  occurs  when  the  step  is 
halfway  between  the  ends  of  the  resonator.  We  accordingly  choose  this 
case.  Furthermore,  we  assume  circular  cross  sectional  areas.  Equation 
(II. A. 11)  then  reduces  to 

tan(^)  -  ^  ,  (III. A. 3) 

2  c  a, 

where  we  have  used  the  dispersion  law  oi  =  ck,  and  where  d,  and  dj  are  the 
diameters  corresponding  to  the  areas  A,  and  Aj,  respectively.  Slapping  the 
other  end  causes  the  boundary  conditions  to  be  switched,  which  inverts  the 
right  side  of  (III. A. 3)  and  thus  leads  to  a  different  frequency.  We  refer 
to  the  two  frecfuencies  as  w.  and  w.,  where  u).  >  w.,  and  to  the  diameters  as 
d.  and  d.,  where  d^  >  d. .  We  wish  to  determine  values  of  d,  and  d.  such 
that  w./o).  is  a  desired  value  in  Table  III.A.l. 

First,  we  must  consider  the  effect  of  the  end  correction.  The  length 
of  the  resonators  is  six  inches,  the  outer  diameter  is  three  inches,  and 
the  characteristic  inner  diameter  is  one  inch.  Because  the  wavelength  is 
much  larger  than  the  wall  thickness  of  the  resonator,  the  end  correction 
amounts  to  the  addition  of  the  length  5L  =  0.305d  to  the  geometrical 
length  L,  where  d  is  the  open-end  diameter  (Pierce,  1989).  We  now  show 
that  the  end  correction  causes  two  errors,  each  of  which  reduces  the  ratio 
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of  the  two  frequencies.  In  (III. A. 3)  the  frequency  is  inversely 
proportional  to  the  length  L.  However,  due  to  the  different  diameters, 
the  values  of  the  effective  length  are  different  for  the  two  cases. 
Because  the  open  end  in  the  case  of  the  frequency  has  a  greater 
diameter  than  that  in  the  u.  case,  the  ratio  will  be  reduced  due  to 
the  end  effect.  This  is  the  dominant  effect.  The  secondary  effect  is  due 
to  the  fact  that  the  end  correction  causes  the  step  to  no  longer  lie 
halfway  between  the  ends  of  the  resonator.  (It  is  nearer  to  the  closed 
end.)  This  also  causes  a  reduction  of  u./o).  because  the  ratio  is  greatest 
when  the  step  is  halfway  between  the  ends.  Due  to  the  fact  that  this  is 
an  extremum,  however,  the  value  of  is  relatively  insensitive  to  the 
location  of  the  step  as  long  as  it  is  near  the  center. 

We  will  include  the  dominant  effect  of  the  end  correction,  but  will 
neglect  the  secondary  effect.  Hence  (III. A. 3)  can  still  be  used  if  the 
length  L  is  replaced  by  L  +  0.305d2.  The  ratio  of  the  frequencies  is  then 
given  by 


w. 


,  0.305d 

‘  *  — T— 
- - 

1  +  - - - 1 


tan-’(.^) 

- ~-J- 

tan  ’  ( ^  ) 


(III. A. 4) 


This  cannot  be  inverted.  However,  because  the  end  effect  is  small,  we  can 
solve  for  values  of  d,  and  d.  by  first  determining  these  values  when  the 
end  correction  is  omitted,  and  then  using  these  as  starting  values  in  a 
numerical  search.  Without  the  end  correction,  (III. A. 4)  is 
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(III. A. 5) 


w, 

oT 


tan-’{-^) 
- T~ 


We  can  invert  this  by  using  the  identity  tan’(l/x)  =  7r/2  -  tan’(x). 
Solving  (III. A. 5)  for  d./d.  then  gives 

-  tan(^_|Lj  ")  ,  (III. A. 6) 

where  A  =  It  should  be  noted  that  the  analytical  determination  of 
d,/d.  without  the  end  correction  is  only  possible  if  the  step  is  located 
halfway  between  the  ends.  otherwise,  a  numerical  solution  is  required. 

Using  (III. A. 6),  we  determine  d./d.  for  the  values  of  A  corresponding 
to  the  musical  intervals  1-3,  1-5,  and  1-8  in  Table  III.A.l.  Next,  we 
uniquely  determine  d^  and  d.  by  arbitrarily  choosing  their  mean  value  to 
be  unity,  which  will  be  taken  to  as  one  inch  in  the  construction  of  the 
resonators.  Finally,  these  values  are  employed  as  a  starting  point  in  a 
simple  numerical  search  to  determine  the  values  of  d,  and  d.  =  2-d^  to  the 
nearest  0.001  that  satisfy  the  more  accurate  equation  (III. A. 4).  This  is 
quickly  accomplished  with  only  a  hand  calculator. 

The  resultant  values  of  d.  and  d.  are  tabulated  in  Table  III. A. 2.  To 
ensure  that  these  diameters  yield  values  of  A  that  are  accurate  to  less 
than  1%,  the  tolerance  of  the  diameters  should  be  ±  0.002  inches.  From 
Table  III.A.l,  a  minor  second  (or  half  step)  corresponds  to  about  a  6% 
increase  in  A,  so  we  consider  deviations  less  than  1%  to  be  acceptable. 

One  may  wonder  if  the  Rayleigh  approximate  result  (II.C.12)  would 
yield  sufficiently  accurate  diameters.  The  best  way  to  investigate  this 
is  to  assume  the  values  of  the  diameters  in  Table  III. A. 2,  substitute  them 
into  the  Rayleigh  result,  determine  the  corresponding  values  of  A  =  w^/w., 
and  compare  these  to  the  musical  values  in  Table  III.A.l.  When  this  is 
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done,  the  deviation  is  found  to  be  0-1%  for  the  1-3  resonator,  0.7%  for 
the  1-5,  and  3.4%  for  the  1-8.  The  Rayleigh  approximation  is  thus 
sufficient  except  in  the  final  case.  It  should  be  noted  that  the  Rayleigh 
approximation  would  be  much  less  accurate  if  we  did  not  employ  the  average 
area  A  =  (A,+A2)/2  as  the  reference  area  in  the  derivation  in  Section  II. C. 
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Table  III.A.l  Relative  frequencies  of  the  notes  of  the  chromatic 
scale. 


Musical  Note 

Relative 

Musical  Note 

Relative 

Frequency 

Frequency 

1 

1 

5 

1.4983 

m 

1.0595 

6 1 

1.5874 

2 

1.1225 

6 

1.6818 

31) 

1.1892 

71) 

1.7818 

3 

1.2599 

7 

1.8877 

4 

1.3348 

8 

2 

4jt 

1.4142 

Table  III. A. 2  Diameters  of  2-pitch  and  3-pitch  resonators  whose  lengths 
are  L  =  6.000  inches. 


Resonator 

Diameter  (inches) 

Diameter  (inches) 

1-3  and  1-2-3 

1.094 

0.906 

1-5  and  1-3-5 

1.164 

0.836 

1-8  and  1-5-8 

1.278 

0.722 
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B.  DESIGN  OF  THE  3-PITCH  RESONATORS 

If  we  imagine  cutting  a  2-pitch  resonator  in  half  along  the  axis,  an 
additional  configuration  (and  therefore  an  additional  frequency)  is  now 
possible.  Furthermore,  slapping  is  no  longer  necessary.  If  one  end  of 
the  resonator  is  set  on  a  smooth  horizontal  surface,  and  the  thumbs  are 
used  as  hinges  where  the  edge  of  the  halves  meet,  we  can  clap  the  halves 
together  to  produce  sound.  These  resonators  are  referred  to  as 
"clappers."  How  can  we  design  these  reaon,<tors  to  approximately  yield 
musical  frequency  intervals? 

The  answer  is  simple.  If  a  slapper  is  imagined  to  be  cut  in  half,  the 
new  frequency  will  approximately  equal  the  average  of  the  original 
frequencies.  Now,  the  musical  notes  2,  3,  and  5  lie  approximately  halfway 
between  1  and  3,  1  and  5,  and  1  and  8,  respectively. 

To  quantitatively  investigate  the  deviation  of  the  middle  frequency 
from  the  desired  musical  value,  we  calculate  the  frequency  ratio 


01 


,  0.305d 

^  *  —z— 

■;-7Tj.305<!; 

L 


TT 


(III.B.l) 


where  d,„,  is  the  root-mean-square  diameter  of  the  end  when  the  diameters 
are  mixed: 


d 


rms 


7 


(I1I.B.2) 


This  is  the  appropriate  value  to  use  in  the  end  correction,  because  it  is 
the  diameter  corresponding  to  the  circle  whose  area  equals  that  of  the 
slapper.  For  L  =  6.000  inches,  and  for  the  slapper  diameters  in  Table 
III. A. 2,  the  deviations  of  the  middle  frequencies  from  the  desired  musical 
values  are  tabulated  in  Table  III.B.l.  Remarkably,  the  deviations  are  all 
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less  than  1%,  which  meets  our  design  criterion.  Hence,  we  will  employ 
this  design.  Of  course,  we  should  not  physically  cut  the  slapper  in  half, 
because  some  material  will  be  loss.  The  effective  diameters  would  then  be 
less  than  the  original  values,  and  the  frequencies  would  shift.  The 
actual  construction  of  the  3-pitch  clappers  is  discussed  in  Section  III.E. 


Table  IlI.B.l  Comparison  of  the  theoretical  middle  frequencies  of  the 
3-Pitch  resonators  to  the  desired  musical  values. 


Resonator 

Relative  Middle 

Desired  Musical 

Deviation 

Frequency  Wo/w, 

Value 

1-2-3 

1.1302 

1.1225 

+  0.69% 

1-3-5 

1.2502 

1.2599 

-  0.77% 

00 

1 

in 

1 

1.5043 

1.4983 

+  0.40% 
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c. 


DESIGN  OF  THE  4-PITCB  RESONATORS 


By  constructing  dissimilar  halves  of  a  clapped  resonator,  we  can 
produce  four  different  frequencies.  It  is  then  possible  co  design  two 
resonators  such  that  one  gives  the  first  four  notes  of  a  musical  scale, 
and  the  other  gives  the  latter  four  notes.  In  this  case,  we  know  from  our 
experience  with  the  previous  resonators  that  the  step  sizes  will  be  small. 
Using  the  Rayleigh  theory  as  a  guide,  we  then  find  that  we  do  not  have  the 
freedom  to  approximate  a  major  scale.  In  fact,  the  symmetry  of  the 
frequency  splittings  (Fig. II. C. 2)  restricts  the  possibilities  among 
standard  scales  to  the  harmonic  minor  (Fig. Ill .C. 1 ) .  In  this  section,  we 
comment  on  a  design  of  a  pair  of  four-pitch  resonators  that  approximately 
yield  this  scale. 

One  way  to  make  the  halves  of  a  resonator  different  vs  construct  the 
steps  to  lie  at  different  locations.  This  would  significantly  increase 
the  difficulty  of  the  theoretical  calculations,  however.  A  simpler  way  is 
to  keep  both  steps  at  the  middle  of  the  resonator,  but  to  alter  the  cross 
sectional  areas  ( Fig. Ill .0. 2 ) .  Equation  (III. A. 4)  thus  still  applies.  As 
encountered  in  Section  III.B,  however,  d.  and  d.  are  now  the  root-mean- 
square  (rms)  values  of  the  semicircular  dxameters  a^  and  b^.  For  the 
configuration  shown  in  Fig. III. C. 2,  d,^  =  (a,^+b,^)/2  and  d.*  =  (a.*+b.*)/2. 
For  the  first  resonator  we  will  choose  diameters  such  that  the  outer  two 
frequencies  differ  by  a  perfect  fourth  (P4)  musical  interval,  and  the 
inner  two  by  a  minor  second  (m2).  For  the  second  resonator,  we  choose  P4 
and  m3  (minor  third).  These  are  diagrammed  in  Fig. III. C. 3.  The  intervals 
between  the  first  and  second  frequencies,  and  the  third  and  fourth,  will 
then  only  approximately  equal  the  desired  musical  values.  We  will  later 
calculate  the  deviations,  and  show  that  they  are  acceptable.  Also,  we 
will  later  calculate  the  lengths  of  the  resonators  such  that  the  second 
resonator  yields  the  latter  four  notes  of  the  scale. 
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To  obtain  the  P4  interval,  we  set 


3 


a. 


(III.C.l) 


where  are  diameters  that  yield  a  P4  interval.  To  obtain  the  m2  (or  m3) 
interval,  we  set 

.2  -  ^  ,  (III.C.2) 


where  are  diameters  that  yield  a  m2  (or  m3)  interval.  Equations 
(III.C.l)  and  (III.C.2)  are  four  equations  in  the  four  unknowns  a,  and  b,. 
However,  they  are  not  consistent  unless  an  area  constraint  is  satisfied: 

*  p  2  -  yi  *  ^2  ,  (III.D.3) 


which  states  that  the  total  area  (of  both  ends)  must  be  independent  of  the 
configuration.  This  is  not  satisfied  if  the  mean  values  of  and  p.,  or 
V,  and  are  chosen  to  be  unity,  as  we  did  before.  We  can  satisfy 

(III.C.3)  by  now  choosing  the  rms  values  to  be  unity: 

•  (III.C.4) 

Equations  (III.C.2)  and  (III.C.3)  then  become  redundant.  In  particular, 
we  are  free  to  conveniently  choose 

-  1  -  •  (III.C.5) 

rms  rms  *  ' 


The  solution  is  then  given  by 


41 


a/  -  „  2  _  1 

bi  -  -  ^  2  ^  I 


(III.C.6) 


By  the  procedure  established  in  Section  III. A,  we  find  =  1.2315  for  a 
P4  interval  and  =  1.1405  for  a  m3,  both  for  a  resonator  of  length  L  = 
6  inches.  We  show  below  that  the  length  of  the  second  resonator  must  be 
L  »  9.142  inches.  We  then  find  =  1.2290  for  a  P4  interval,  and  = 
1.0465  for  a  m2.  The  resultant  values  of  the  diameters  a^  and  b^  for  both 
resonators  are  tabulated  in  Table  III.C.l. 

We  have  now  designed  two  4-pitch  resonators  whose  outer  and  inner 
frequency  intervals  are  correct.  We  must  determine  how  accurate  are  the 
intervals  between  the  first  and  second  frequencies,  as  well  as  the  third 
and  fourth.  This  can  be  accomplished  as  in  Section  III.B.  The  results 
are  that  the  frequencies  deviate  by  about  1%  from  the  musical  values.  As 
stated  in  Section  III.B,  we  consider  this  to  be  acceptable  because  the 
notes  of  a  m2  interval  (the  smallest)  deviate  by  6%. 

Finally,  we  must  determine  lengths  of  the  two  resonators  such  that  the 
second  will  yield  the  notes  5, 6k, 7,  and  8  relative  to  the  first,  which  has 
the  notes  1,2,3k,  and  4  (Fig. III. C. 3)  Hence,  we  demand  that 

-  0f,  ,  (III.C.7) 


where  0  =  1.4983  (which  corresponds  to  a  perfect  fifth).  The  frequency 
splittings  of  each  resonator  are  symmetrical  about  the  central  values  f, 
and  f|, .  These  are  the  frequencies  of  the  unperturbed  (uniform  cross 
sectional  area)  resonators.  Because  the  outer  frequencies  of  each 
resonator  have  the  Seune  ratio,  it  must  be  that  f^/f»  =  fj/f,  or 
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(III.C.8) 


^6 


f// 


T, 


Now,  the  ratio  of  the  unperturbed  frequencies  is 


f„  L,  *  SL 
7,  ~  L„  *  SL  ' 


(III.C.9) 


where  L,  and  L,,  are  the  lengths  of  the  resonators,  and  6L  is  the  end 
correction  (0.305  multiplied  by  the  open-end  diameter).  Combining 
(III.C.7),  (III.C.8),  and  (III.C.9),  we  find  that  we  must  choose  the 

lengths  of  the  resonators  such  that 

0-1^  • 


We  set  L„  =  6.000  inches  and  approximate  the  diameter  as  1.000  inch.  The 
end  correction  is  thus  6L  =  0.305.  Solving  for  L,  gives 

If,  —  ~  1  /  6lf  ,  (III. C. 11) 

which  yields  L,  =  9.142  inches.  The  effect  of  the  end  correction  here  is 
about  2%,  which  is  not  negligible  according  to  our  1%  criterion. 
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log  (frequency) 


_  8  _  8  _  8 

_  7  _  7 

_  7b 

_  6 

_  6b  _  6b 

_  5  _  5  _  5 

_ 4  _  4  _  4 

_  3 

_  3b  _  3b 

_  2  _  2  _  2 

_  1  _  1  _  1 

Major  Natural  Harmonic 

Minor  Minor 


Fig.III.C.l  Standard  musical  scales.  The  4-pitch  resonator  is 
conducive  to  the  harmonic  minor  scale. 
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Fig. III. C. 2  Geometry  of  a  4-pitch  resonator.  The  quantities  a,  and 
b,  are  diameters. 
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log  frequency 


Resonator  II 


1  ^ 

—  8 

y 

1 

m3 

1 

1 

6b 

II 
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r 

1  6 

_ 

V 

P4 


Resonator  I 


3b 

2 


A  m2 

V 


P4 


I 

V 


Pig. III. C. 3  Frequency  diagram  for  the  2-pitch  resonators  that  yield 
a  harmonic  minor  scale. 
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Table  III.C.l  Dimensions  of  a  set  of  4-pitch  resonators  that  yield 
a  harmonic  minor  scale. 


Resonator 

Length 
( inches ) 

Diameters 
( inches) 

Diameters 
( inches ) 

l-2-3l)-4 

9.142 

1.129,  0.851 

1.087,  0.904 

5-6l)-7-8 

6.000 

1.171,  0.792 

1.045,  0.953 

D .  MATERIAL 

The  desire  to  obtain  experimental  results  that  would  best  demonstrate 
the  theory,  as  previously  discussed  in  Chapter  II,  established  a  criterion 
for  the  basic  requirements  of  this  thesis.  To  achieve  our  criterion, 
acoustic  resonators  had  to  be  constructed  out  of  a  material  which  would 
produce  resonators  of  sufficient  quality  to  meet  our  needs.  There  were 
several  factors  in  the  selection  of  the  construction  material  for  the 
resonators.  The  most  important  factor  was  the  desire  to  produce  the 
cleanest  signal  and  clearest  tone  possible  without  the  use  of  an  excessive 
excitation  force.  The  use  of  an  excessive  excitation  force  could  lead  to 
over  ringing  and  reverberation  of  harmonics  of  the  resonant  frequency. 
The  result  of  an  excessive  excitation  force  can  be  seen  in  Fig.III.D.l. 
The  beating  of  the  resonant  frequency  is  clearly  evident.  Additionally, 
the  cost,  availability,  and  inherent  manufacturing  difficulties  of  the 
different  materials  (as  will  be  discussed  in  Section  III.E),  were  also 
contributing  factors  in  our  decision. 

Several  acoustic  resonators  of  both  the  two-pitch  (slapper)  and  three- 
pitch  (clapper)  variety  were  constructed  out  of  different  materials  in 
order  to  determine  the  best  material  for  this  research.  The  two  different 
methods  of  exciting  the  resonators  (Section  I.B)  provided  a  means  to 
determine  which  of  the  materials  would  be  acceptable  for  use  in 
constructing  the  resonators.  Based  on  tests  conducted  with  slappers 
constructed  from  acrylic  and  wood  (fir),  the  tone  quality  was  found  not  to 
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be  dependent  on  the  material  of  the  slapper,  as  expected.  Additionally, 
the  quality  of  the  free  decay  time  series  plot  were  also  similar. 
However,  the  fir  proved  to  be  a  more  difficult  material  to  machine  than 
the  acrylic.  This  was  apparent  when  a  comparison  of  the  fir  slapper  to 
the  acrylic  slapper  showed  that  the  critical  interior  walls  of  the  fir 
resonator  did  not  bore  as  smoothly.  Since  the  tolerances  of  the  interior 
diameters  were  ±  0.002  inches  (Section  III. A),  sanding  of  the  interior 
wall  smoother  was  not  a  viable  option.  Furthermore,  from  a  demonstration 
viewpoint,  the  clear  acrylic  allows  one  to  view  the  step  in  the  resonator. 
Based  on  these  considerations,  the  material  of  choice  for  the  two-pitch 
slappers  was  acrylic. 

The  excitation  signal  of  the  clappers  was  created  from  the  hitting  (or 
clapping)  of  the  two  symmetric  halves  flush  onto  each  other  (Section  I.B). 
Here,  unlike  the  slappers,  the  resulting  excitation  signal  is  dependent  on 
the  material  of  the  clapper.  The  first  clapper  produced  was  machined  from 
a  stock  of  aluminum.  The  higher  density  of  the  aluminum  allowed  the 
excitation  of  usable  signals  with  only  a  minimum  use  of  force.  However, 
the  tone  quality  of  the  resonator  was  poor  due  to  the  highly  metallic  ring 
of  the  signal.  As  with  the  slappers,  the  second  choice  for  a  construction 
material  was  acrylic.  The  acrylic  proved  to  be  a  very  acceptable  choice 
of  material  for  construction.  This  was  due  to  the  fact  that  the 
excitation  of  a  clean  signal  with  purer  tone  quality,  using  a  minimum 
amount  of  force,  was  more  easily  reproduceo. 

The  need  to  insure  that  a  more  acceptable  material  was  not  overlooked 
resulted  in  a  clapper  constructed  from  maple.  While  the  maple,  being 
denser  than  fir,  made  for  easier  machining,  the  maple  clapper  did  not 
produce  a  signal  as  clean  nor  with  as  long  of  a  useable  record  length  as 
the  acrylic  clappers.  Comparison  of  Fig. III. D. 2  (a  1-2-3  resonator  made 
from  maple)  and  Fig. III. D. 3  (a  1-2-3  resonator  constructed  from  acrylic) 
clearly  demonstrates  the  better  signal  quality  of  acrylic  versus  maple. 
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The  tone  quality  of  the  maple  clapper  also  was  inferior  to  the  acrylic. 
These  factors  again  pointed  toward  the  acrylic  as  the  material  of  choice. 
It  should  be  noted  that  several  more  exotic  woods  such  as  rosewood  and 
mahogany  were  considered  but  due  to  the  exorbitant  cost  of  these  exotic 
woods  they  were  disregarded.  These  would  probably  yield  a  sound  similar 
to  the  maple.  Therefore,  acrylic  was  selected  to  be  the  construction 
material  for  all  of  the  resonators. 
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Fig. III. 0.1  Example  of  beating  of  signal  due  to  use  of  an  excessive 
excitation  force  for  the  1-8  slapper  resonator. 
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Fig. III. D. 2  Time  series  plot  of  a  12-3  clapper  resonator  constructed 
from  maple.  Note  the  deterioration  of  the  signal  after  40  ms. 
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Fig. III. 0.3  Time  series  plot  of  the  1-2-3  clapper  constructed  from 
acrylic.  Note  that  the  acrylic  produces  a  usable  signal  for  the 
entire  record  length. 
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E.  CONSTRUCTION  OP  RESONATORS 

A  high  degree  of  quality  and  precision  was  required  in  the 
construction  of  the  resonators.  In  order  for  the  experimental  results  not 
to  deviate  more  than  ±  1.0%  from  the  desired  musical  value,  a  tolerance  of 
±  0.002  inches  on  the  interior  diameters  of  the  resonators  was  stipulated. 
To  achieve  this  tolerance  the  following  method  was  used  to  construct  the 
resonators. 

Initially,  the  two,  three,  and  four-pitch  resonators  were  started 
using  the  same  process,  utilizing  a  blank  acrylic  solid  cylinder,  seven 
inches  long,  cut  from  a  three  inch  diaimeter  acrylic  rod.  The  blank  was 
then  placed  on  a  three  chuck  lathe  and  a  pilot  hole,  whose  diameter  was 
less  than  the  smallest  diameter  (required  for  that  particular  resonator) , 
was  drilled  through  the  cylinder.  To  obtain  the  desired  tolerance,  the 
requisite  diameters  were  bored  out  using  a  boring  bit  attached  to  a  boring 
bar.  It  should  be  noted  that  the  diameter  of  the  pilot  hole  would  have  to 
be  large  enough  to  accept  the  boring  bar.  Such  that  it  was  critical  that 
the  bored  cavities  be  concentric,  during  the  boring  phase  the  cylindrical 
blank  was  not  removed  from  the  lathe  and  both  diameters  were  bored  from 
the  larger  diameter  side.  In  the  construction  of  ail  of  our  resonators, 
the  step  change  from  one  diameter  to  the  other  occurred  at  the  halfway 
point  of  the  resonator.  This  completed  the  construction  for  the  two-pitch 
resonators. 

The  three-pitch  resonators  were  actually  two  separate  but  identical 
two-pitch  resonators  cut  in  half  lengthwise.  Each  two-pitch  resonator  was 
cut  longitudinally  off  centered  and  then  milled  down  to  the  desired 
exterior  diameter,  such  that  the  two  halves  were  concentric.  The  interior 
diameters  of  the  two  and  three-pitch  resonators  are  given  in  Table 
III. A. 2. 
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The  four-pitch  resonators  were  constructed  in  the  same  manner  as  the 
three-pitch  resonators,  except  that  instead  of  two  identical  halves  each 


two-pitch  resonator  was  of  different  interior  diameters.  The  dimensions 
for  the  four-pitch  resonator  are  given  in  Table  III.C.l. 
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IV.  EXPERIMENT 


In  this  chapter,  we  discuss  our  experiment  to  test  the  theory.  In  the 
first  two  sections,  we  discuss  the  methods  of  detecting  and  exciting  the 
sound.  In  the  last  two  sections,  we  present  the  data  and  its  analysis. 

A.  DETECTION 

In  this  section,  the  methods  used  to  obtain  data  will  be  discussed. 
This  will  include  the  equipment  used  and  factors  which  could  lead  to 
experimental  or  random  errors.  Our  desire  was  to  capture  a  high  quality 
free-decay  signal  for  each  of  the  frequencies  of  every  resonator.  Since 
each  excitation  of  a  resonator  is  an  independent  event,  we  required  the 
ability  to  store  multiple  signals  of  a  single  resonator,  for  a  detailed 
analysis  at  a  later  time.  Additionally,  without  the  use  of  an  anechoic 
chamber,  we  wanted  to  ensure  that  the  captured  free-decay  signal  of  each 
resonator  was  a  true  representation  of  that  free-decay  without  the  effects 
of  outside  influences.  The  detailed  analysis  included  measurement  of  the 
period  for  each  individual  cycle,  a  computer  generated  FFT  of  for  every 
decay  signal,  and  a  determination  of  the  quality  facto-  for  each  decaying 
signal . 

There  were  several  stipulated  considerations  involved  in  the  obtaining 
of  measurements  for  this  thesis,  that  provided  the  criterion  for  our 
method.  The  basic  requirement  of  determining  the  resonant  frequencies 
for  each  of  two,  three,  and  four-pitch  resonators  to  within  ±  0.2  Hz  was 
the  primary  driving  factor.  To  capture  these  signals,  we  would  require 
just  the  basics  of  acoustic  signal  measuring.  Since  the  resonant 
frequency  was  our  desired  goal,  we  needed  some  type  of  an  oscilloscope 
capable  of  displaying  the  decaying  signal  in  such  a  manner  that  an 
accurate  measurement  of  the  periods  could  be  made.  Further  equipment 
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required  included  a  combination  high  pass/low  pass  filter  to  filter  out 
undesired  signals,  and  a  microphone  sensitive  enough  to  respond  to  the 
excitation  signal.  However,  the  microphone  must  not  be  too  sensitive, 
because  we  desire  that  the  influence  of  background  noise  be  kept  to  a 
minimum. 

The  first  piece  of  equipment  selected  was  the  oscilloscope.  The 
requirement  to  have  the  capability  of  storing  numerous  records  made  the 
Nicolet  Digital  Oscilloscope  (NIC-310)  the  ideal  choice  for  our  work.  The 
NIC-310  provided  us  the  capability  to  record  one  separate  signal  after 
another  on  a  3.5"  floppy  disk  and  then  at  a  later  time  retrieve  the  signal 
for  analysis.  The  digital  oscilloscope  would  sample  4000  data  points  for 
each  single  experimental  run,  regardless  of  the  length  of  signal  to  be 
captured.  Additionally,  the  NIC-310  has  the  capability  to  expand  the  x 
and  y  axis  of  a  stored  signal,  on  the  screen,  thus  allowing  a  precise 
measurement  of  each  of  the  4000  data  points,  from  the  digital  numerical 
display.  After  some  initial  experimentation  using  the  NIC-310  while 
attempting  to  record  data  from  a  1-5  slapper  resonator  (an  original 
version  without  end  effect  corrections  included),  the  average  length  of  a 
good  quality  time  series  of  the  data  signal  was  approximately  80  ms. 
Thus,  tor  the  measurement  of  all  resonators  the  time  per  point  setting 
used  on  the  NIC-310  was  20  ^e.  Of  particular  note  is  that,  since  the 
method  used  to  measure  the  periods  was  to  mark  the  zero  crossings  of  the 
X  axis  for  each  cycle,  the  DC  offset  on  the  NIC-310  needed  to  be  zero.  It 
was  discovered  that  the  NIC-310  had  a  DC  voltage  drift  of  approximately  1 
mV  per  minute.  Therefore,  the  NIC-310  required  a  zeroing  of  the  DC  offset 
voltage  after  each  experimental  run.  The  selection  of  a  acceptable 
microphone  proved  to  be  more  complicated  than  the  selection  of  the 
oscilloscope.  The  first  microphone  used  was  a  Realistic  #33-918  Dynamic 
Microphone.  In  the  early  stages  of  our  work  the  Realistic  appeared  to  be 
more  than  adequate  for  our  measurements.  However,  later  it  was  discovered 
that  mass  element  of  the  Realistic  produces  a  ring  down  effect.  This 
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effect  produces  a  frequency  of  approximately  200  to  300  Hz,  with  higher 
harmonics  in  the  600  and  1500  Hz  range  measured.  This  revelation 
suggested  that  use  of  the  Realistic  microphone  could  cause  serious 
problems  with  the  reliability  of  the  recorded  data.  Since,  the  ring  down 
frequency  (fundamental  and  harmonics)  falls  in  the  middle  of  the 
theoretical  frequency  range  of  all  the  resonators  of  300  to  700  Hz,  we 
therefore  disregarded  the  Realistic  microphone.  Our  next  selection  was 
the  General  Radio  microphone  #1962-9610  with  a  #1560-P42  preamplifier  and 
a  #1560-P62  power  supply.  Test  revealed  that  since  the  General  Radio  is 
of  the  electret  type,  the  microphone  did  not  produce  the  undesired  ring 
down  effect.  However,  the  General  Radio  power  supply  did  pick  up  and 
amplify  60  Hz  noise.  This  60  Hz  noise  occurs  even  though  the  power  supply 
was  on  batteries,  if  the  electrical  power  cord  was  not  connected  to  the 
120  V  wall  outlet.  This  was  attributed  to  the  lack  of  proper  60  Hz 
grounding  on  the  building  occupied  by  the  lab.  Therefore,  in  the  course 
of  this  experiment  the  power  supply  was  left  plugged  into  the  wall  outlet 
while  the  power  supply  was  running  on  batteries,  to  ground  out  the  60  Hz 
interference.  Two  additional  types  of  microphones  were  tested  to  insure 
that  the  General  Radio  was  indeed  the  best  for  our  purposes.  A  Bruel- 
Kjaer  type  was  tested  but  proved  to  be  too  sensitive  to  background  noises 
to  make  it  usable  in  our  lab.  Finally,  a  microphone  was  constructed  using 
a  Realistic  Electret  element,  which  cost  all  of  one  dollar.  Upon 
comparison  of  test  with  the  Bruel-Kjaer  and  the  General  Radio  microphones, 
the  Electret  produced  results  that  match  up  fairly  well.  However,  since 
numerous  measurement  were  to  be  carried  out  over  an  extended  period,  we 
decided  that  the  General  Radio  microphone  was  better  suited  to  our 
purposes.  Since,  the  General  Radio  microphone  provided  us  with  a  higher 
degree  of  reliability. 

The  final  piece  of  equipment  selected  was  the  high  pass/low  pass 
filter.  Two  different  filters  were  considered  the  Ithaco  1201  and  the 
Krhon-Hite  3202R.  Since,  the  General  Radio  preamplifier  has  the  ability 
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to  produce  a  gain  of  lOX  (the  desired  gain  to  produce  a  reasonable 
amplitude  to  be  measured  on  the  NIC-310)  the  built  in  gain  on  the  Ithaco 
would  not  be  required.  We  tested  both  filters  for  internal  noise  and 
cutoff  frequencies  with  the  low  pass  filter  set  at  3000  Hz  and  the  high 
pass  set  at  150  Hz.  It  was  concluded  that  either  filter  would  have  been 
adequate  for  our  purposes,  but  our  decision  to  use  the  Krohn-Hite  filter 
was  based  on  the  sharper  cutoff  frequencies.  A  diagram  of  the 
experimental  setup  is  shown  in  Fig. IV. A. 1. 

The  experimental  method  used  to  obtain  the  data  was  first  to  excite 
the  resonator  repetitively  until  a  desirable  signal  was  displayed  on  the 
oscilloscope,  then  save  the  signal  on  a  3.5”  floppy  disk.  After  a  number 
of  signals  consider  to  be  a  fair  representation  of  the  decay  signal  were 
acquired,  the  best  three  individual  signals  were  then  recalled  to  the 
screen,  the  x  and  y  axis  were  expanded  and  the  measurement  of  period  for 
each  cycle  was  made  to  the  nearest  0.01  ms  (using  the  method  of  zero 
crossings  as  discussed  previously).  The  time  series  plot  of  a  typical 
decay  signal  (in  this  case  the  1-3  slapper)  is  shown  by  Fig. IV. A. 2.  After 
the  periods  for  each  cycle  were  recorded,  the  average  period  over  the 
entire  length  of  the  signal  was  calculated.  From  this  average  period  the 
average  resonant  frequency  was  determined.  The  individual  period  values 
plotted  in  sequential  order  from  time  zero  for  one  data  run  of  the  1-3 
slapper  can  be  seen  in  Fig. IV. A. 3.  The  record  data  stored  on  a  disk  was 
then  used  to  produce  a  computer  generated  FFT  of  the  time  series  using 
PCMATLAB  as  seen  in  Fig. IV. A. 4  (for  the  same  record  as  before).  The 
frequency  resolution  of  12.5  Hz  for  the  FFT  was  predetermined  from  the 
sampling  frequency  of  the  NIC-310.  The  sampling  frequency  used  for  our 
work  was  50  kHz.  This  rather  high  frequency  was  necessary,  due  to  the 
short  time  span  of  our  signal,  in  order  to  obtain  sufficient  data  points 
to  reconstruct  the  signals  and  measure  the  periods  within  the  desired 
accuracy.  Our  desire  to  have  a  computer  generated  FFT  of  the  decay  signal 
was  to  observed  the  frequency  shifting  of  the  overtones  due  to  the 
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nonuniformity  of  the  resonators.  The  arrows  in  Fig. IV. A. 4  represent  the 
overtones  for  the  uniform  resonator  (i.e.  a  resonator  with  no  step). 
Lastly,  because  it  is  not  clear  that  the  decay  should  be  linear  (Section 
V.B),  we  examined  the  logarithm  of  the  voltage  verses  time.  The  decays 
were  indeed  found  to  be  linear,  and  the  quality  factor  Q  for  each 
resonator  was  determined  using 
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(IV. A. 1) 


where  oi  is  the  resonant  frequency  and  a  is  the  slope  of  the  natural 
logarithm  of  the  absolute  value  for  the  voltage  amplitudes  plotted  verses 
time.  Fig.  IV.  A.  5  shows  the  plot  used  to  determined  a,  for  a  sample 
resonator. 
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Fiq.lV.A.2  Typical  time  series  of  the  decay  signal  for  a  2-pitch 
"slapper"  resonator.  This  time  series  plot  is  the  high  frequency  for 
a  1--3  resonator. 


Fig. IV. A. 3  Dispersion  plot  of  the  individual  periods  for  each  cycle 
of  a  typical  2-pitch  "slapper”  resonator.  This  plot  is  of  the  high 
frequency  for  a  1-3  resonator. 
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Fig. IV. A. 4  Computer  generated  FFT  plot  using  PCMATLAB  for  a  typical 
2-pitch  "slapper"  resonator.  Note  the  frequency  shifts  of  the 
overtones  of  the  nonuniform  case  as  compared  to  the  overtones  for  the 
uniform  case  (arrows). 
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Fig. IV. A. 5  Plot  of  the  logarithm  of  the  absolute  voltage  amplitude 
verses  time.  This  plot  was  used  to  compute  the  Quality  factor  for  a 
typical  2-pitch  "slapper"  resonator. 


B .  EXCITATION 


It  should  be  noted  that  for  demonstration  purposes,  as  discussed  in 
Section  I.B,  that  the  palm  of  the  hand  can  be  used  for  excitation  of  the 
two-pitch  resonator.  However,  in  reality  the  hand  excitation  of  the 
slapper,  while  producing  usable  data,  proved  to  be  a  difficult  method  to 
reproduce  continuously  and  required  numerous  experimental  runs  to  acquire 
just  one  good  signal.  This  was  discovered  to  be  a  result  of  a  changing  in 
the  boundary  conditions  at  the  end  of  the  slapper  caused  by  the  palm  of 
the  hand  not  sealing  the  resonator  completely.  Additionally,  after  data 
comparison,  it  was  observed  that  a  small  fraction  of  the  palm  of  the  hand 
entered  into  the  interior  of  the  resonator,  thereby  changing  the  acoustic 
length  of  the  slapper.  This  changing  of  the  acoustic  length  of  the 
slapper  was  unacceptable  for  our  measurements,  since  this  produced  a 
change  in  the  resonant  frequency  that  would  not  be  reproducible  from  data 
run  to  data  run.  Therefore,  a  different  method  of  exciting  the  slapper 
was  required  that  produced  the  same  boundary  conditions  every  time.  It 
was  observed  that  if  a  #12  rubber  stopper  was  used  to  excite  the  slapper 
by  striking  the  stopper  against  the  end  of  the  slapper,  the  output  signal 
was  not  only  easier  to  reproduce  time  and  time  again  but  that  the  signal 
was  exceptionally  clean.  This  clean  signal  meant  that  the  entire  length 
of  the  signal  from  time  zero  to  80  ms  was  now  usable.  By  comparison  of 
Fig.IV.B.l,  a  time  series  of  the  1-8  slapper  using  excitation  by  the  palm 
of  the  hand  and  Fig. IV. B. 2  where  excitation  was  by  use  of  a  #12  stopper, 
it  can  be  seen  that  the  signal  in  Fig.  IV.  B.  2  is  more  desirable  for 
experimental  purposes. 

For  the  case  of  excitation  of  the  three  and  four-pitch  resonators  the 
method  described  in  Section  I.B  was  effective  in  producing  good  quality 
data . 
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Fig.IV.B.l  Times  series  decay  plot  of  the  low  frequency  for  a  1-8 
slapper  resonator  excited  by  using  the  palm  of  the  hand. 


66 


d 

•  d 

d 

—  O 

d 

(yi(OA)  3?eifO/\ 

5 

C'J  ^ 

=  9  ? 

Fig. IV. B. 2 

Times 

series 

plot  of  the 

low 

frequency  for  a  1-8  slapper 

resonator  excited  using  a  #12  rubber  stopper. 


67 


c. 


DATA 


In  this  section  the  data  acquired  and  reduced  as  discussed  in  Section 
IV. A  will  be  presented.  For  each  resonator,  three  separate  experimental 
runs  were  used  for  the  data  reduction  of  each  particular  frequency.  After 
the  initial  decay  signal  was  recorded  the  methods  discussed  in  Section 
IV. A  were  used  to  obtain  the  results  included  in  Table  IV.C.l  through 
Table  IV. C. 8.  To  clarify  the  presented  data,  the  average  periods  were 
measured  using  an  average  of  40  consecutive  cycles  in  the  case  of  the 
higher  frequencies  of  the  resonators  and  an  average  of  28  consecutive 
cycles  for  the  lower  frequencies.  From  the  average  period  the  average 
resonant  frequency  was  then  calculated.  The  quality  factor  was  then 
determined  in  the  manner  discussed  in  Section  IV. A. 

As  can  be  observed  from  Table  IV.C.l  through  Table  IV. C. 8,  there  is 
relatively  good  agreement  between  the  average  resonant  frecfuencies  for 
each  experimental  run  for  a  given  resonator.  Nominally,  the  difference  in 
resonant  frequencies  appears  to  less  than  0.4  Hz,  this  is  as  expected  due 
to  the  experimental  errors  encountered,  which  would  affect  each 
independent  run  to  some  degree.  The  factor  that  proved  to  be  the  most 
sensitive  and  caused  large  deviations  in  the  resonant  frequencies  between 
each  experimental  run  was  the  air  temperature.  The  temperature  had  little 
effect  when  the  accuracy  was  required  to  be  within  one  hertz.  However, 
when  accuracy  was  required  to  be  v'ithin  tenths  of  hertz  the  temperature 
was  the  limiting  factor.  To  overcome  this  problem  required  the  cycling 
through  the  frequencies  of  each  resonator  by  capturing  one  record  for  a 
particular  frequency  then  capturing  one  record  for  the  next  frequency. 
This  method  was  repeated  until  three  records  had  been  captured  for  each 
frequency. 
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This  method  proved  to  be  successful  and  allowed  the  average  resonant 
frequencies  to  be  computed  from  experimental  runs  whose  deviations  were  in 
the  tenths  of  hertz.  The  quality  factors  for  the  all  of  the  resonators 
were  approximately  in  the  40  to  60  range. 
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Table  XV.C.l  Data  from  the  1-3  Slapper  Resonator 


Record 

Period  Average 
(ms) 

Resonant  Frequency 
Average  (Hz) 

Quality 

Factor 

S1345 

2.136 

468.3 

51 

S1346 

2.136 

468.2 

59 

S1348 

2.135 

468.5 

52 

S3141 

1.698 

588.9 

52 

S3142 

1.698 

589.0 

54 

S3144 

1.697 

589.2 

50 

Table  IV. C. 2  Data  from  the  1-5  Slapper  Resonator 


Record 

Period  Average 
(ms) 

Resonant  Frequency 
Average  (Hz) 

Quality 

Factor 

S1525 

2.385 

419.2 

44 

S1529 

2.386 

419.2 

47 

S1530 

2.386 

419.1 

43 

S5131 

1.594 

627.3 

66 

S5133 

1.593 

627.7 

69 

S5135 

1.594 

627.6 

64 

Table  IV. C. 3  Data  from  the  1-8  Slapper  Resonator 


Record  Period  Average  Resonant  Frequency  Quality 
(ms) Average  (Hz) Factor 


51822  _ 2.895 _ 345.4 _ ^ 

51823  _ 2.895 _ 345.4 _ « 

S1803  2.894  345.5  40 


58104  _ 1.453 _ 688.3 _ ^ 

58105  _ 1.454 _ 688.5 _ 48 

58106  1.453  688.4  47 


Table  IV. C. 4  Data  from  the  1-2-3  Clapper  Resonator.  The  records 
correspond  to  the  lower,  middle  and  high  frequency,  respectively. 


Record 

Period  Average 
(ms) 

Resonant  Frequency 
Average  (Hz) 

Quality 

Factor 

C12353 

2.158 

463.4 

41 

C12355 

2.155 

464.0 

42 

C12356 

2.157 

463.6 

43 

C23157 

1.910 

523.6 

45 

C23159 

1.910 

523.5 

43 

C23160 

1.912 

523.2 

44 

C32146 

1.719 

581.7 

40 

C32147 

1.718 

581.9 

40 

C32162 

1.718 

581.8 

39 

Table  ZV.C.S  Data  from  the  1-3-5  Clapper  Resonator.  The  records 
correspond  to  the  lower,  middle  and  high  frequency,  respectively. 


Record 

Period  Average 
(ms) 

Resonant  Frequency 
Average  (Hz) 

Quality 

Factor 

C13526 

2.391 

418.3 

50 

C13527 

2.391 

418.3 

50 

C13529 

2.392 

418.1 

48 

C35131 

1.925 

519.4 

41 

C35132 

1.925 

519.6 

42 

C35133 

1.926 

519.2 

42 

C53136 

1.601 

624.6 

39 

C53149 

1.602 

624.4 

36 

C53150 

1.602 

624.4 

36 
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Table  IV. C. 6  Data  from  the  1-5-8  Clapper  Resonator.  The  records 
correspond  to  the  lower,  middle  and  high  frequency,  respectively. 


Record 

Period  Average 
(ms) 

Resonant  Frequency 
Average  (Hz) 

Quality 

Factor 

C15872 

2.892 

345.8 

44 

C1S875 

2.890 

346.0 

41 

C15878 

2.891 

345.9 

45 

C58173 

1.963 

509.4 

53 

C58176 

1.965 

508.9 

52 

C58179 

1.963 

509.4 

50 

C85174 

1.454 

687.8 

39 

C85177 

1.453 

688.2 

34 

C85180 

1.454 

687.8 

40 
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Table  IV. C. 7  Data  from  the  l-2-3b-4  Clapper  Resonator.  The  records 
correspond  from  the  lowest  to  the  highest  frequency,  respectively. 


Record 

Period  Average 
(ms) 

Resonant 
Frequency 
Average  (Hz) 

Quality  Factor 

C4120 

3.311 

302.0 

36 

C4124 

3.310 

302.1 

36 

C4128 

3.308 

302.3 

41 

C4221 

2,929 

341.4 

42 

C4225 

2.932 

341.1 

45 

C4229 

2.931 

341.2 

42 

C4322 

2.759 

362.4 

38 

C4326 

2.761 

362.2 

36 

C4330 

2.761 

362.2 

38 

C4423 

2.484 

402.6 

38 

C4427 

2.486 

402.3 

38 

C4431 

2.485 

402.4 

39 

t 
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Table  IV. C. 8  Data  from  the  5-6 k -7-8  Clapper  Resonator.  The  records 
correspond  from  the  lowest  to  the  highest  frequency,  respectively. 


Record 

Period  Average 
(ms) 

Resonant 
Frequency 
Average  (Hz) 

Quality  Factor 

C4505 

2.226 

449.2 

49 

C4509 

2.225 

449.4 

43 

C4513 

2.226 

449.2 

56 

C4602 

2.090 

478.5 

42 

C4610 

2.089 

478.7 

43 

C4614 

2.090 

478.5 

41 

C4715 

1.754 

570.1 

52 

C4716 

1.755 

569.8 

54 

C4717 

1.755 

569.8 

57 

C4804 

1.671 

598.4 

47 

C4808 

1.671 

598.4 

44 

C4818 

1.671 

598.4 

48 
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D.  ANALYSIS 

In  this  final  section  of  the  experimental  chapter  the  final  results 
from  the  research  conducted  are  displayed  in  Table  IV.D.l  through  Table 
IV. D. 3.  The  average  frequencies  are  calculated  by  taking  the  average  of 
the  three  resonant  frequencies  (Table  IV.C.l  through  Table  IV. C, 8) 
measured  for  each  particular  frequency  of  each  resonator.  The  relative  ■< 

frequency  is  simply  the  ratio  of  the  high  frequency  to  the  low  frequency 
which  is  the  musical  interval  or  step  of  the  resonator. 

It  can  clearly  be  seen  that,  of  the  resonators  constructed  for  testing 
by  using  the  design  criteria  given  by  Section  III. A  through  Section  III.C 
for  the  two,  three,  and  four-pitch  resonators,  respectively,  the  relative 
frequencies  were  in  good  agreement  with  the  desired  musical  intervals.  It 
should  also  be  noted  that  by  incorporating  end  effects  into  our  design  of 
the  resonators,  we  were  able  to  get  the  experimental  data  to  agree  with 
the  desired  musical  value  to  within  our  predetermined  criterion  of  less 
than  roughly  one  percent  deviation. 
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Table  IV.D.l  Final  Results  for  the  2-Pitch  Resonators 


Resonator 

Low  Frequency 
(Hz) 

High 

Frequency 

(Hz) 

Relative 

Frequency 

Deviation 
from  Musical 
Value 

1-3 

468.3 

589.0 

1.258 

-  0.2  % 

1-5 

419.2 

627.5 

1.497 

-  0.1  % 

1-8 

345.4 

688.4 

1.993 

-  0.4  % 

Table  IV. D. 2  Final  results  for  the  3-Pitch  Resonators 


Resonator 

Low 

Frequency 

(Hz) 

Higher 

Frequency 

(Hz) 

Relative 

Frequency 

Deviation 
from  Musical 
Value 

1-2-3 

463.7 

523.4 

1.129 

+  0.6  % 

581.8 

1.255 

-  0.4  % 

1-3-5 

418.2 

519.4 

1.242 

-  1.4  % 

624.5 

1.493 

-  0.4  % 

1-5-8 

345.9 

509.2 

1.473 

-  1.7  % 

687.9 

1.989 

-  0.6  % 
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V.  CONCLUSIONS  AND  FUTURE  WORK 


In  this  chapter,  we  summarize  the  results  and  then  discuss  possible 
future  research  involving  nonuniform  resonators. 

A.  CONCLUSIONS 

We  have  investigated  piecewise  uniform  acoustic  resonators  whose 
configurations  can  be  simply  and  quickly  altered,  leading  to  changes  in 
the  frequency.  Various  resonators  have  been  constructed  such  that  the 
frequencies  correspond  to  musical  notes.  There  are  two  types  of 
resonators:  single-piece  “slappers" ,  whose  ends  are  slapped  in  order  to 
generate  the  sound,  and  double-piece  "clappers",  which  are  clapped 
together  while  one  set  of  ends  rests  on  a  smooth  surface.  The  apparatus 
are  primarily  useful  as  educational  demonstrations  of  the  effects  of 
nonuniformity  on  standing  wave  modes.  The  demonstrations  have  the 
advantage  of  being  inexpensive  and  simple  to  perform. 

On  an  elementary  level,  the  demonstrations  can  be  appreciated  as  an 
effect  of  symmetry-breaking.  The  explanation  is  not  elementary,  however. 
We  have  presented  three  theoretical  approaches;  matching  solutions  of  the 
wave  equation,  applying  the  energy  method  of  Rayleigh,  and  applying  the 
theorem  of  adiabatic  invariance.  The  second  and  third  approaches,  which 
are  perturbative  and  yield  the  same  result,  offer  physical  explanations. 
The  first  approach,  which  is  more  accurate,  was  employed  in  the  design  of 
the  resonators.  We  obtained  excellent  agreement  with  this  theory  when  end 
effects  were  approximately  included,  and  observed  small  but  significant 
deviations  from  the  perturbative  theory. 


B.  FUTURE  WORK 

There  are  several  possible  avenues  for  future  research  regarding 
piecewise  uniform  resonators.  In  this  section,  we  briefly  discuss  some  of 
these. 

The  standard  theory  is  not  exact,  although  our  data  show  no  indication 
of  a  breakdown  of  the  theory.  It  would  be  interesting  to  perform  a  t 

careful  experiment  to  attempt  to  observe  deviations.  It  may  be  that  the 
breakdown  only  occurs  when  the  length  of  the  resonator  becomes  comparable 
to  the  diameter;  that  is,  when  the  system  is  no  longer  approximately  one¬ 
dimensional.  Analytically,  one  could  begin  with  the  exact  two-dimensional 
transmission  step  theory  of  Morse  and  Ingard  (1968).  This  may  be 
relevant . 

We  are  not  aware  of  any  quantitative  experiment  in  the  breakdown  of 
adiabatic  invariance.  A  piecewise  uniform  resonator  with  a  movable  step 
or  length  may  be  a  feasible  candidate.  As  discussed  in  Section  II. F,  when 
the  discontinuity  in  area  is  large,  pairs  of  modal  freqxiencies  can  be 
separated  by  a  small  aunount  (Fig. II . F.  1) .  Exciting  one  mode  and  then 
altering  the  geometry  such  that  the  frequencies  of  this  and  another  mode 
become  close,  may  lead  to  a  partial  excitation  of  the  nearby  mode.  This 
would  be  a  breakdown  of  adiabatic  invariance.  It  may  inherently  be 
connected  to  a  breakdown  of  the  constancy  of  E/w  for  the  original  mode. 

Finally,  we  consider  a  question  regarding  the  free  decay  of  a  mode  of 
a  closed-closed  piecewise  uniform  resonator.  Due  to  the  different  areas, 
the  cpaality  factors  of  the  individual  uniform  sections  will  differ.  One 
may  then  wonder  whether  or  not  a  mode  decays  as  a  pure  exponential.  * 

Specifically,  consider  a  large  widening  step  halfway  along  the  resonator 

f 

(Aj  »  A,  and  a  =  1/2  in  Fig. II. A. 1).  By  the  results  of  Section  II. F,  the 
fundamental  and  first  overtone  have  nearby  frequencies.  Furthermore,  the 
results  suggest  that  the  behavior  can  be  represented  as  a  system  of  two 
coupled  identical  oscillators  with  different  damping  parameters.  It  is 
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not  difficult  to  show  that,  according  to  this  model,  the  decay  of  either 
normal  mode  is  not  purely  exponential.  However,  the  effect  is  significant 
only  if  the  coupling  is  weak  (or  if  the  frequencies  of  the  modes  deviate 
by  a  small  amount).  It  would  be  interesting  to  analytically  and 
experimentally  investigate  the  free  decay  of  a  sound  mode  in  the  geometry 
described  above.  Our  experience  is  that  the  typical  intuition  of  an 
acoustician  is  that  the  decay  is  always  purely  exponential,  regardless  of 
the  geometry. 
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